Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 482,982 A→G 19.9% N33D (AAT→GAT)  USA300HOU_0461 → possible PAP2 superfamily phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730482,9820AG19.9% 37.2 / 4.9 20N33D (AAT→GAT) USA300HOU_0461possible PAP2 superfamily phosphatase
Reads supporting (aligned to +/- strand):  ref base A (7/9);  new base G (2/2);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

ACATTTGTAAAATGAAAATGATGATACGAGGAGATTGTAAGATGATAGATAAAAAATTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAATTACTATATGACGTTTAGTGCGGAAATTGGAGATGTC  >  CP000730/482845‑483116
                                                                                                                                         |                                                                                                                                      
aCATTTGTAAAAGGAAAATGATGATACCAGAAAATTCAAAAAAGAAAAAAAAAAAATTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATc                                                                                                                                     <  2:88000/141‑1 (MQ=255)
                 atgatgATACGAGGAGATTGTAAGATGATAGATAAAAAATTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAat                                                                                                                    >  2:55672/1‑141 (MQ=255)
                     tgatACGAGGAGATTGTAAGATGATAGATAAAAAATTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTCACAAATCAAGAATGGCTTAAAAATATAg                                                                                                                <  1:158519/141‑1 (MQ=255)
                                   tGTAAGATGATAGATAAAAAATTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTa                                                                                                  >  2:188905/1‑141 (MQ=255)
                                               gataAAAAATTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACa                                                                                      >  2:241300/1‑141 (MQ=255)
                                                  aaaaaaTTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGat                                                                                   <  1:213266/141‑1 (MQ=255)
                                                   aaaaaTTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCTCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAGATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGatt                                                                                  >  2:136662/1‑141 (MQ=255)
                                                   aaaaaTTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAGATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGatt                                                                                  >  2:144567/1‑141 (MQ=255)
                                                                      aaTGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGTGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCacgacg                                                               >  2:168350/1‑141 (MQ=255)
                                                                       aTGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGACTATTTCGGTGAGCCACAACGt                                                              <  1:128098/141‑1 (MQ=255)
                                                                                  tttATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTaa                                                   <  1:188603/141‑1 (MQ=255)
                                                                                  tttATTGTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTaa                                                   <  2:159908/141‑1 (MQ=255)
                                                                                               gcgcTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAAt                                      <  1:188905/141‑1 (MQ=255)
                                                                                                   tGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAGATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAATTACt                                  <  1:144567/141‑1 (MQ=255)
                                                                                                   tGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAGATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTACGGTGAGCCACAACGTCAGTATGTTAACATTTTGTTTAATTACt                                  <  1:136662/141‑1 (MQ=255)
                                                                                                                  gTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAATTACTATATGACGTTTAGTg                   <  1:55672/141‑1 (MQ=255)
                                                                                                                          tACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAATTACTATATGACGTTTATTGCGGAAAtt           >  2:90611/1‑141 (MQ=255)
                                                                                                                          tACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAATTACTATATGACGTTTAGTTCGGAAAtg           >  2:73989/1‑140 (MQ=255)
                                                                                                                          tACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGATTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAATTACTATATGCCGTTTAGTGCGGCAAtt           >  2:47014/1‑141 (MQ=255)
                                                                                                                                   gttaCAAATCAAGAATAGCATAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAATTACTATATGACGTTTAGTGCGGAAATTGGAGATGTc  <  1:94051/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
ACATTTGTAAAATGAAAATGATGATACGAGGAGATTGTAAGATGATAGATAAAAAATTAACATCACCGAAAATGACAGTGCCTTTATTTTTAATCGCGCTGATTGTATTTATAGGTATGTTTTACAGTGTAGTGACAAATCAAGAATGGCTTAAAAATATAGATATGGGATCATTAACATGGTTTACAGATTATTTCGGTGAGCCACAACGTCAGTATGTTAACAATTTGTTTAATTACTATATGACGTTTAGTGCGGAAATTGGAGATGTC  >  CP000730/482845‑483116

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