Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,428 T→C 16.5% K82E (AAA→GAA)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4280TC16.5% 95.9 / 13.0 48K82E (AAA→GAA) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (24/16);  new base C (6/2);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

TTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGT  >  USA300TCH1516_ALE/1744288‑1744563
                                                                                                                                              |                                                                                                                                       
ttCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAAttt                                                                                                                                         <  1:139159/141‑1 (MQ=255)
    aTCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCct                                                                                                                                     >  2:165982/1‑141 (MQ=255)
         acctccAGTACCTAATAAATCATCTGTAATTAACACACTTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGTTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCttt                                                                                                                                <  2:358661/141‑1 (MQ=255)
          ccaccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt                                                                                                                               <  1:171755/141‑1 (MQ=255)
             ccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTa                                                                                                                            >  1:402665/1‑141 (MQ=255)
                    cTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTg                                                                                                                     >  1:126303/1‑141 (MQ=255)
                    cTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCACTAAACCGACATTGCCCCCATCCCAGCGCAGACGCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTATTCTACCAGATg                                                                                                                     >  2:301548/1‑141 (MQ=255)
                                cTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATcc                                                                                                         <  2:171885/141‑1 (MQ=255)
                                cTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTATAACAGGTGCTAAGCCAATcc                                                                                                         >  1:354602/1‑141 (MQ=255)
                                cTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATcc                                                                                                         >  1:228883/1‑141 (MQ=255)
                                             acGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGc                                                                                            <  2:402665/141‑1 (MQ=255)
                                               gTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACAAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCCAAGCCAATCCCCCTTCAAAAATCtc                                                                                          >  2:38225/1‑140 (MQ=255)
                                                  gACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAg                                                                                       >  2:160951/1‑141 (MQ=255)
                                                    ccTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGa                                                                                     <  1:244817/141‑1 (MQ=255)
                                                     ttGGTGTAATTGCATCGTTGTGGATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAGGAc                                                                                    <  2:114456/140‑1 (MQ=255)
                                                       ggTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAGGACAg                                                                                  <  2:15882/141‑1 (MQ=255)
                                                           tAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACATGACAGccca                                                                              >  1:392332/1‑139 (MQ=255)
                                                             aTTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAAGGTCATACGGCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATg                                                                            <  1:199019/141‑1 (MQ=255)
                                                             aTTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATTAGCTACAGGACAGCCAATg                                                                            >  1:233695/1‑141 (MQ=255)
                                                             aTTGCATCTTTGTGCATTGTTAAAAACTATCAAACCTATGCGAGGTCATA‑‑CTTATAACGAATGACTTCACGAGGTAATTCACCTTCTTTTCTAACACATTCAAAGCCAAACACCATTAAATAATCTACAGTACAGACaaaa                                                                            >  2:92179/1‑139 (MQ=255)
                                                              ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGa                                                                           >  1:106707/1‑141 (MQ=255)
                                                              ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGa                                                                           >  1:203517/1‑141 (MQ=255)
                                                                 cATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATCATaa                                                                        >  2:204245/1‑141 (MQ=255)
                                                                 cATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaa                                                                        >  1:184581/1‑141 (MQ=255)
                                                                   tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGCTGCTAAGCAAATCCCCATTGAATAAGCTACAGGACAGCCAACGATaaaa                                                                      >  2:317883/1‑140 (MQ=255)
                                                                        gtgCATTGGTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAAGGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg                                                                 <  2:411987/141‑1 (MQ=255)
                                                                         tgCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcgc                                                                <  1:236876/141‑1 (MQ=255)
                                                                             ttcttAAAACATTTGTAACATATTCTAGGTCATC‑‑CTCATAACGAATGACTTCACGAGGTACTTTCTCTTCTTTCCTAACAGGTGCAAAGCCAATCCCCATTGAATACGCTGCGGCACATCCAATGATAAAGCCACGCGCtt                                                            <  2:321808/138‑1 (MQ=255)
                                                                                     aCATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCt                                                    >  1:354982/1‑141 (MQ=255)
                                                                                     aCATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCATGCGCTTCAGGTCCt                                                    >  2:383061/1‑141 (MQ=255)
                                                                                      cATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCACGTCCtc                                                   >  1:343927/1‑140 (MQ=255)
                                                                                       aTTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTAc                                                  <  1:165982/141‑1 (MQ=255)
                                                                                         ttGTACCATGTTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACaa                                                <  2:355746/141‑1 (MQ=255)
                                                                                           gTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAAGTCCTACAACg                                              >  1:102692/1‑141 (MQ=255)
                                                                                              ccATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAACCTACAGGACAGCCAATGATAAAGCCACGCCCCTCAGGTCCTACCACGcta                                           >  2:314545/1‑141 (MQ=255)
                                                                                                         gTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGTTCCTACAACGATATCAACATctct                                >  2:25126/1‑141 (MQ=255)
                                                                                                          tCATA‑‑CTCATAACGAGTGACTTCACGCGGTAATTTCCCTTCTTTTCTAACGGGTGCTAAGCCAATCCCCTTTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTg                               <  1:166933/141‑1 (MQ=255)
                                                                                                           cATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCTAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGt                              >  1:408337/1‑141 (MQ=255)
                                                                                                                  tCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACCATATCAACATCCCTGTCtttt                         >  2:337259/1‑141 (MQ=255)
                                                                                                                   cATAACGAATGACTTCACGCGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtgttt                        >  1:417548/1‑140 (MQ=255)
                                                                                                                           aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACCGCCAATGATAAAGCCACGCGCTTCAAGTCCTACAACGATATCAACATCTCTGTCTTTTACGTATTCt                >  2:55499/1‑141 (MQ=255)
                                                                                                                           aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                <  1:21604/141‑1 (MQ=255)
                                                                                                                           aTGACTTCACGAGGTAATTTCCCTTCTATTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                >  1:229502/1‑141 (MQ=255)
                                                                                                                            tGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCCTTTGCGTATTCTa               >  1:74056/1‑141 (MQ=255)
                                                                                                                              aCTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACa             >  1:154813/1‑141 (MQ=255)
                                                                                                                                ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAt           <  1:381501/141‑1 (MQ=255)
                                                                                                                                       ggTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCt    >  1:65142/1‑141 (MQ=255)
                                                                                                                                         tAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt  <  2:126303/141‑1 (MQ=255)
                                                                                                                                         tAATCTACCTTGTTTTCTAACGAGTGCAATGGTAATCCCCTGTTATTAAGCTACAGGACAGCTAATTATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt  <  2:6529/141‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                       
TTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATA‑‑CTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGT  >  USA300TCH1516_ALE/1744288‑1744563

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: