Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 665,392 T→G 100% L79* (TTA→TGA)  paiA_1 → Spermidine/spermine N(1)‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE665,3920TG97.2% 68.6 / ‑6.0 21L79* (TTA→TGA) paiA_1Spermidine/spermine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (14/6);  minor base C (1/0);  total (15/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACAACAGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATGAAGAA  >  USA300TCH1516_ALE/665255‑665528
                                                                                                                                         |                                                                                                                                        
gACAACAGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATgata                                                                                                                                       >  1:637616/1‑141 (MQ=255)
gACAACAGTTAATCACTTGTTAGAAGCTTAATATAATGAACATCATTTATAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAGTgata                                                                                                                                       >  1:157202/1‑141 (MQ=255)
 acaacaGTTCATCACTTGTTAGAAGCTTCATCTCTGGTGCAGCTGGTCAAGAAAAGACTTCAAGGACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGatat                                                                                                                                      >  1:596676/1‑141 (MQ=255)
  caacaGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATAtt                                                                                                                                     >  1:502704/1‑141 (MQ=255)
              aCTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCAtt                                                                                                                         <  1:282485/141‑1 (MQ=255)
              aCTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCAtt                                                                                                                         <  1:770591/141‑1 (MQ=255)
               cTTGTTAGAAGCTTCGTATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTAAttt                                                                                                                        >  1:603413/1‑140 (MQ=255)
                   ttAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGtt                                                                                                                    >  1:177992/1‑141 (MQ=255)
                   ttAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACGTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGtt                                                                                                                    >  1:133129/1‑141 (MQ=255)
                        aGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAcc                                                                                                               <  1:366057/141‑1 (MQ=255)
                                          tcatTTAAAGAAAAGAATTCAAGAACAATTATTCTTAGGCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAgg                                                                                             >  1:455146/1‑141 (MQ=255)
                                                   gaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTAc                                                                                    <  1:807400/141‑1 (MQ=255)
                                                              aaGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCCTTGTTaa                                                                         >  1:403543/1‑141 (MQ=255)
                                                                    aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTACCTTTATTTACGGTGACGAAACTTTCTTATCAGCTAATTATGTTAATCCAGAATCGCAACATACAGTTTATGGTACAGCATTGTTAAATGAAg                                                                   >  1:801615/1‑141 (MQ=255)
                                                                     attagtCGTAGTCGGTGAAGAAGTTAATGACATCGTTGGCTTTGCTAACTTTATTTACCGTGAAGAATGAGATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAGTGTTAAATGAAgg                                                                  <  1:480291/141‑1 (MQ=255)
                                                                               gTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGt                                                        >  1:424115/1‑141 (MQ=255)
                                                                                         aagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATa                                              >  1:145291/1‑141 (MQ=255)
                                                                                                 gACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGaa                                      >  1:716548/1‑141 (MQ=255)
                                                                                                                        tatttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAGTTAACGCTGACAAGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGataa               <  1:354994/141‑1 (MQ=255)
                                                                                                                          tttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGataata             >  1:7040/1‑141 (MQ=255)
                                                                                                                             aCGGCGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCCCCACACACCGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATaaaa          >  1:555803/1‑141 (MQ=255)
                                                                                                                                   aagaaTGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGGTAATAAAAATgaag    >  1:314954/1‑141 (MQ=255)
                                                                                                                                     gaaTGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATgaagaa  >  1:721107/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GACAACAGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATGAAGAA  >  USA300TCH1516_ALE/665255‑665528

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