Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,856,087	T→A	P80P (CCT→CCA)	USA300HOU_2698 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,856,087	0	T	A	96.4%	84.9 / ‑2.7	28	P80P (CCT→CCA)	USA300HOU_2698	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (17/10);  total (18/10)
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.68e-01

GAGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAATTCA  >  CP000730/2855950‑2856224                                                                                                                                          |                                                                                                                                          gagaTGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCaaaa                                                                                                                                        <  1:73477/141‑1 (MQ=255)   gaTGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTACAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACg                                                                                                                                      >  1:283250/1‑141 (MQ=255)     tGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTc                                                                                                                                    >  2:271103/1‑141 (MQ=255)           gCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATTCGCCAAAACGTCAACaaa                                                                                                                              >  2:147795/1‑141 (MQ=255)               tGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGAt                                                                                                                          >  1:268210/1‑141 (MQ=255)                 gTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTAAATATACAGTTAACAGGCGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAg                                                                                                                        >  1:331719/1‑141 (MQ=255)                   ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                      <  2:243426/141‑1 (MQ=255)                      aaaGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAtttt                                                                                                                   >  2:18848/1‑141 (MQ=255)                           tataTCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgt                                                                                                              <  1:315165/141‑1 (MQ=255)                               tccaatattgatataGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgttgtt                                                                                                          >  1:53940/15‑141 (MQ=255)                               tCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgttgtt                                                                                                          <  1:159475/141‑1 (MQ=255)                                        gAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAATCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTATATAGTTGTTACAgcagca                                                                                                 >  2:80680/1‑141 (MQ=255)                                         aaCAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAt                                                                                                <  2:221133/141‑1 (MQ=255)                                                     ttCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAat                                                                                    >  1:194448/1‑141 (MQ=255)                                                                       ttAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTa                                                                  >  1:66593/1‑141 (MQ=255)                                                                               gAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGaa                                                          <  2:243847/141‑1 (MQ=255)                                                                                         catcatTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGt                                                >  2:2888/1‑141 (MQ=255)                                                                                                 ccTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCt                                        >  2:88181/1‑141 (MQ=255)                                                                                                   ttttAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGCTAGATTTAATAGTTGATACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTgg                                      <  1:312609/141‑1 (MQ=255)                                                                                                           tAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTg                              >  2:137937/1‑141 (MQ=255)                                                                                                            aaCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGa                             >  2:73886/1‑141 (MQ=255)                                                                                                            aaCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGa                             <  1:297020/141‑1 (MQ=255)                                                                                                               aGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACat                          >  1:312177/1‑141 (MQ=255)                                                                                                                       tgtCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGAc                  >  1:302029/1‑141 (MQ=255)                                                                                                                         tCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCt                <  2:331719/141‑1 (MQ=255)                                                                                                                         tCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTATATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGATTTGACATATACGGACCt                >  2:352364/1‑141 (MQ=255)                                                                                                                               aaTCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAata          <  1:169165/141‑1 (MQ=255)                                                                                                                                  cATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTa       >  2:47378/1‑141 (MQ=255)                                                                                                                                       gcCAAAACGTCAACAAATGATAGATTTTATAGTTGGTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATAGACGGACCTCAAATATTAATTCa  <  2:285665/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                          GAGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAATTCA  >  CP000730/2855950‑2856224

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