Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,616,078 G→T 10.5% P113Q (CCA→CAA)  nudF ← ADP‑ribose pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,616,0780GT10.5% 88.7 / 2.4 38P113Q (CCA→CAA) nudFADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (19/15);  new base T (2/2);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

TTCATTAGCATCGATTTAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAGCAATGGTTTTT  >  USA300TCH1516_ALE/1615946‑1616216
                                                                                                                                    |                                                                                                                                          
ttCATTAGCATCGATTTAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACca                                                                                                                                    >  1:147977/1‑141 (MQ=255)
   aTTAGCATCGATTTAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACcatac                                                                                                                                 <  2:202828/141‑1 (MQ=255)
       gCATCGATTTAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACatat                                                                                                                             <  1:174422/141‑1 (MQ=255)
             aTTTAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTacaa                                                                                                                       >  1:31591/1‑141 (MQ=255)
             aTTTAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTCTTCATCGCAAAAACCTGGTGAACCATACATATCTacaa                                                                                                                       >  2:285881/1‑141 (MQ=255)
              tttAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGCACCATACATATCTacagc                                                                                                                      >  2:63373/1‑139 (MQ=255)
              tttAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCACCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTacaac                                                                                                                      <  2:20197/141‑1 (MQ=255)
               ttAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTacaaca                                                                                                                     >  1:247260/1‑141 (MQ=255)
                    aTTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATgtgt                                                                                                                <  2:249907/141‑1 (MQ=255)
                     ttCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGACCTGTACCTTTTGCTAAATTATCCGGAAAATATATTGATAATTGTTCATCGCAACCACCTGGTGAACCATACACTTCTACACCATGCGtt                                                                                                               >  2:367866/1‑141 (MQ=255)
                       ctcAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTaa                                                                                                             >  1:259601/1‑141 (MQ=255)
                         cAATCGGAACTTTAATGCCCAAGACGAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAAct                                                                                                           <  2:205047/141‑1 (MQ=255)
                                aaCTTTAATGACTTCGACAAATTTATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGc                                                                                                    <  2:204930/141‑1 (MQ=255)
                                         gACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATcc                                                                                           <  1:148011/141‑1 (MQ=255)
                                               gACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATAGCCTGTttc                                                                                     >  1:4712/1‑141 (MQ=255)
                                                 cAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACCTGTGTTAACTCTTTGGCAATATATCCTGTttctt                                                                                   >  2:253627/1‑141 (MQ=255)
                                                  aaaTTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTTGCAATATATCCTGTttcttc                                                                                  >  2:391336/1‑141 (MQ=255)
                                                     ttcatcttcatcTAAGTGAACTGTACCTTCTTCTAAATTATCCGTACAATCTCTTAATCATTGTTCATCGCAAACACCTGGTGAACCATACATATCTACCACACGTGTTAACTCTTTAGCCATATATCCTGTttcttttcc                                                                               >  2:280526/1‑137 (MQ=255)
                                                                      gAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCg                                                              <  1:270067/141‑1 (MQ=255)
                                                                             aCCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAAACATGCATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGAc                                                       >  1:290588/1‑141 (MQ=255)
                                                                               cttcttCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTTGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACTTTTTGCCGCTTCAActt                                                     >  2:169316/1‑140 (MQ=255)
                                                                               cttcttCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTTGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGAttc                                                     >  2:168331/1‑141 (MQ=255)
                                                                                     ctAAATTATCCGTAAAATATATTGATAATTGTTCATCGCCAAAACCTGATGAACCATACCTCTCTCCAACCTGTGTTCACTCTTTGGCGATATATCCTGTTTCTTCTTCTAATTACCGTTTTGCCGCTTCTACTCTATCtt                                               >  1:131161/1‑141 (MQ=255)
                                                                                      tAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCTGCTTCGACTCTATCTtc                                              >  1:116855/1‑141 (MQ=255)
                                                                                      tAAATTATCCGTAAAATATATTGACACTTTTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTtc                                              <  2:315901/141‑1 (MQ=255)
                                                                                        aaTTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTtcat                                            >  1:177380/1‑141 (MQ=255)
                                                                                                             aTAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACTACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTg                       <  2:292192/141‑1 (MQ=255)
                                                                                                                   gTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAAttt                 >  1:135951/1‑141 (MQ=255)
                                                                                                                   gTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATATTCATCATCTTCTAATTTACCAGCTGGAAttt                 >  1:33098/1‑141 (MQ=255)
                                                                                                                     tCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTCAACTCTTTGGCACATTATCCTGTCTCTTCTTCTAATTCCCGTTTTGCCGCTTCGCCTCTATCTTTTTAATCTTCTAATTTACCAACTCCACTTTCt               >  2:115266/1‑141 (MQ=255)
                                                                                                                       aTCGCAAAAACCTTGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAg             <  1:169316/141‑1 (MQ=255)
                                                                                                                       aTCGCAAAAACCTTGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAg             <  1:168331/141‑1 (MQ=255)
                                                                                                                       aTCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAg             <  1:31027/141‑1 (MQ=255)
                                                                                                                            aaaatCCTGGTGAACCATACATATCTACGACATGTGTTAACTCTTTTGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAAGTTCTAGCAATg        <  2:282859/141‑1 (MQ=255)
                                                                                                                             aaaaCCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAGCAATgg       <  1:238989/141‑1 (MQ=255)
                                                                                                                                aCCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAGCAATGGttt    <  2:183979/141‑1 (MQ=255)
                                                                                                                                aCCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAGCAATGGttt    >  1:231467/1‑141 (MQ=255)
                                                                                                                                  cTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAGCAATGGttttt  <  1:63373/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                          
TTCATTAGCATCGATTTAACATTCTCAATCGGAACTTTAATGACTTCGACAAATTCATCTTCATCTAAGTGAACTGTACCTTCTTCTAAATTATCCGTAAAATATATTGATAATTGTTCATCGCAAAAACCTGGTGAACCATACATATCTACAACATGTGTTAACTCTTTGGCAATATATCCTGTTTCTTCTTCTAATTCACGTTTTGCCGCTTCGACTCTATCTTCATCATCTTCTAATTTACCAGCTGGAATTTCTAGCAATGGTTTTT  >  USA300TCH1516_ALE/1615946‑1616216

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