Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,607 T→C 100% D22G (GAT→GGT)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,6070TC88.2% 69.8 / 2.1 33D22G (GAT→GGT) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base C (16/10);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTGAAAAATCTTGATA  >  USA300TCH1516_ALE/1744471‑1744743
                                                                                                                                        |                                                                                                                                        
aGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCtt                                                                                                                                      <  1:441695/141‑1 (MQ=255)
    aCAGGACAGCCAATGATAAAGCCACGCGCATCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATACTTCTAGTAATAACCTTGcaa                                                                                                                                  >  2:283625/1‑141 (MQ=255)
     cAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAAc                                                                                                                                 <  2:244467/141‑1 (MQ=255)
              cAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTg                                                                                                                        >  1:323226/1‑141 (MQ=255)
              cAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATAACAACATCTCTGTCTTTTGCGTCTTCTACAATTTTATCAGTTTCATCGCCATCAGCTTCAACACAATCACTAATTGTCGTAATAAACTTCTAACTAccacgtc                                                                                                                        >  2:128067/1‑138 (MQ=255)
                 tGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCACCATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGtt                                                                                                                     >  1:416334/1‑141 (MQ=255)
                     aaaGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATACTTCTAGTAATACCCTTGACTCTAACATATGGttgatt                                                                                                                 >  1:231009/1‑137 (MQ=255)
                              cgcTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCTTTGAAACTAACACCTGGTTTCGGCCAATCTTg                                                                                                        >  1:246851/1‑141 (MQ=255)
                              cgcTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTg                                                                                                        <  1:41125/141‑1 (MQ=255)
                                 ttCAGGTCCTACAACGATATCAACATCTCTGTCATTTGCGTATTCTACAATTTTATCTGTTGCCTAGCAATATGCTTCACCATTATACATAATTGTAGTAATATCCTTGAAACTAACAACTGGTTTCCGCCAATCTTCAAc                                                                                                     >  2:204873/1‑141 (MQ=255)
                                  tCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACt                                                                                                    <  1:421705/141‑1 (MQ=255)
                                    gggTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATACTCGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTc                                                                                                  <  2:486122/140‑1 (MQ=255)
                                             aaCGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTa                                                                                         <  1:27611/141‑1 (MQ=255)
                                                  tatCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCt                                                                                    >  1:12727/1‑141 (MQ=255)
                                                   atCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTCTTGCtt                                                                                   >  2:217927/1‑141 (MQ=255)
                                                    tCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCttt                                                                                  >  1:166370/1‑141 (MQ=255)
                                                      aaCATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTaa                                                                                >  2:262886/1‑141 (MQ=255)
                                                       ccctctcTTTCATTGGCGTATGTTACCATTTGAGCCGTTGAACAGCCATATTCTTCACTATTATCCATAATTGTAGTAATAACCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACTTAGTGCTTTaaa                                                                               <  2:260404/138‑1 (MQ=255)
                                                                 cTTTTGCTTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACTCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAata                                                                     <  1:38562/141‑1 (MQ=255)
                                                                 cTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAata                                                                     >  1:155796/1‑141 (MQ=255)
                                                                                       aTCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAACCTAACCCCTGGTTTCGGCCAATCTTGAACTTCTGATCCGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTcccca                                               >  2:405201/1‑137 (MQ=255)
                                                                                       aTCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTCATATTTCCTCCTAAATTGCGCACGa                                               >  1:184627/1‑141 (MQ=255)
                                                                                       aTCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGa                                               <  1:197194/141‑1 (MQ=255)
                                                                                          tGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGCACTTTGGAGACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACaa                                            >  1:18598/1‑141 (MQ=255)
                                                                                              gCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATtgt                                        <  1:332740/141‑1 (MQ=255)
                                                                                                aTAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGCATTGCTTTCAATCCATTAATATTTCCTCCTCAATTGCTCACTACAAttatca                                      >  1:481061/1‑137 (MQ=255)
                                                                                                aTAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAAAAACTTTGGTCGGGCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGa                                      >  1:71633/1‑141 (MQ=255)
                                                                                                        atGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATcc                              >  1:185090/1‑141 (MQ=255)
                                                                                                                        cATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTGa              <  1:247839/141‑1 (MQ=255)
                                                                                                                               gtagtaATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCACTTTTTTATTTCTGAAGActct       >  2:184099/1‑141 (MQ=255)
                                                                                                                               gtagtaATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTGAAAAATCt       >  1:259910/1‑141 (MQ=255)
                                                                                                                                 agtaATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATACAATTTTTTATTTCTGAAAAATCttt     >  2:324777/1‑140 (MQ=255)
                                                                                                                                    aaTATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTTTGAAATCTTGata  >  1:311535/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                        
AGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTGAAAAATCTTGATA  >  USA300TCH1516_ALE/1744471‑1744743

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