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
*	CP000730	482,982	0	A	G	19.9%	37.2 / 4.9	20	N33D (AAT→GAT)	USA300HOU_0461	possible 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: ‑