Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	108,148	A→G	11.8%	K1037R (AAA→AGA)	USA300HOU_0106 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	108,148	0	A	G	11.8%	86.9 / 4.7	34	K1037R (AAA→AGA)	USA300HOU_0106	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (9/21);  new base G (2/2);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTGG  >  CP000730/108009‑108288                                                                                                                                            |                                                                                                                                             tATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGaaa                                                                                                                                             <  2:154309/141‑1 (MQ=255)   ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAAcc                                                                                                                                           <  1:56091/141‑1 (MQ=255)   ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAAcc                                                                                                                                           <  1:443731/141‑1 (MQ=255)   ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTGGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATATAGATGAAAcc                                                                                                                                           <  1:60695/141‑1 (MQ=255)       gTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCAttt                                                                                                                                       >  1:432536/1‑141 (MQ=255)       gTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCAttt                                                                                                                                       >  1:293801/1‑141 (MQ=255)                     taataCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAGACCATTTTATGAGACGAtttt                                                                                                                         >  1:71913/1‑141 (MQ=255)                     taataCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAGACCATTTTATGAGACGAtttt                                                                                                                         >  1:302309/1‑141 (MQ=255)                         ccccAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa                                                                                                                     <  2:359129/140‑1 (MQ=255)                         aCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa                                                                                                                     <  2:84766/141‑1 (MQ=255)                         aCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa                                                                                                                     <  2:269728/141‑1 (MQ=255)                           ccAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaaga                                                                                                                   >  2:189597/1‑141 (MQ=255)                              aGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaagaaag                                                                                                                >  1:227057/1‑141 (MQ=255)                                 tGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAgaaagaaa                                                                                                             <  1:143909/141‑1 (MQ=255)                                  ggTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAAt                                                                                                            <  1:442946/141‑1 (MQ=255)                                   gTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATg                                                                                                           <  1:216517/141‑1 (MQ=255)                                   gTGCTATGAAAAGGTCATGAGAAAAACCCAAATTGTTAAAAGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAACGAATACAGCTGAAACCATCTTATGAGACGATTTTTAAAGCAAGAAATg                                                                                                           <  2:298691/141‑1 (MQ=255)                                                gTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATAc                                                                                              <  1:460406/141‑1 (MQ=255)                                                gTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATAc                                                                                              <  2:372951/141‑1 (MQ=255)                                                             ccAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGtatata                                                                                 <  2:432536/141‑1 (MQ=255)                                                             ccAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAAAAAAAAAAAAAAAAGtatata                                                                                 <  2:293801/141‑1 (MQ=255)                                                                 aCTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAgg                                                                             >  1:221888/1‑141 (MQ=255)                                                                  cTTGTTAAACGTTGCAGTGACAAGAGGAAAGAGAGAGGTGTGTGGGATGGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGa                                                                            >  2:20480/1‑141 (MQ=255)                                                                             ttGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAAAAAAAAAAAAAGTATATAGAGGAAGTTATACAtt                                                                 <  2:364321/141‑1 (MQ=255)                                                                             ttGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAAAAACAAAAAAAAAAGTATATAGAGGAAGTTATACAtt                                                                 <  2:418026/141‑1 (MQ=255)                                                                              tGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACAttt                                                                <  1:98540/141‑1 (MQ=255)                                                                                    gACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAg                                                          <  2:221888/141‑1 (MQ=255)                                                                                     aCAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAgg                                                         <  2:317835/141‑1 (MQ=255)                                                                                     aCAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAgg                                                         <  2:348154/141‑1 (MQ=255)                                                                                                                  tGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTaac                            <  1:298101/141‑1 (MQ=255)                                                                                                                      gACATGCAAAGAATACAGATGAGACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACt                        <  2:302309/141‑1 (MQ=255)                                                                                                                      gACATGCAAAGAATACAGATGAGACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAAAAAAAAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACt                        <  2:71913/141‑1 (MQ=255)                                                                                                                            cAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGatat                  >  1:161637/1‑141 (MQ=255)                                                                                                                              aaGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTa                >  2:178362/1‑141 (MQ=255)                                                                                                                                    aCAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCa          >  1:378502/1‑141 (MQ=255)                                                                                                                                            aaCCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTgg  <  1:177476/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                             TATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTGG  >  CP000730/108009‑108288

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