Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	522,999	A→G	14.5%	M1M (ATG→GTG) †	USA300HOU_0485 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	522,999	0	A	G	14.5%	64.4 / 2.8	27	M1M (ATG→GTG) †	USA300HOU_0485	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/9);  new base G (2/2);  total (16/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 = 1.00e+00

ATATATTAGACATTTCTATGCAGGTGAAGCTTGAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAATGATGAACCA  >  CP000730/522864‑523139                                                                                                                                        |                                                                                                                                             atatatTAGACATTTCTATGCAGGTGAAGCTTGAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGt                                                                                                                                         >  1:244117/1‑141 (MQ=255)    aatTACACATTTCTATGCAGGTGAAGCGTGAGGGACATGAACAGCCACGTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAAGAATTATTACATTTGGATCGCCATGAAATATTTGAAAAAATAATGCGTTTa                                                                                                                                      <  2:75612/140‑1 (MQ=255)              ttCTATGCAGGTGAAGCTTGAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAGTGCGTTTAGAaatgaatg                                                                                                                            >  2:202352/1‑141 (MQ=255)              ttCTATGCAGGTGAAGCTTGAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAGTGAGTTTAGAaaggaaag                                                                                                                            >  2:270772/1‑139 (MQ=255)                 tATGCAGGTGAAGCTTGAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAGTGtca                                                                                                                         >  2:152402/1‑141 (MQ=255)                    gCAGGTGAAGCTTGAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGtcaatc                                                                                                                      <  2:164474/141‑1 (MQ=255)                     cAGGTGAAGCTTGAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGtcaatca                                                                                                                     >  1:14364/1‑141 (MQ=255)                              cTTAAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCaaa                                                                                                            <  2:104125/141‑1 (MQ=255)                                  aGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGaa                                                                                                        <  1:115794/141‑1 (MQ=255)                                      aCATGAACAGCCACTTGAATAAAAAATAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACtt                                                                                                    <  1:79971/141‑1 (MQ=255)                                      aCATGAACAGCCAATTGAATATAAATTAGAAGAAATAGAAATTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACtt                                                                                                    <  1:283477/141‑1 (MQ=255)                                             cAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAAAAATGAATGTCAATCAACATTAAAACGAAACTGCAGCAtt                                                                                             >  2:263470/1‑141 (MQ=255)                                                 cACTTGAATATAAATAAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAGAGGAAACTTCAGAATTAAAg                                                                                         >  1:210057/1‑141 (MQ=255)                                                                 agaagaAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAAt                                                                         >  1:142074/1‑141 (MQ=255)                                                                                 aTGCATTCAGGAGGCATTATTACATATGGATCGCAATGAAATATTTGAAAAAATAATGAGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGAAGTTGAATTAGTAGAAGAAAATGt                                                         >  2:140492/1‑141 (MQ=255)                                                                                            aggCATGATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATCAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCATTTGAATTAGTAGAAGAAAATGTAGCGCTTCAAc                                              <  2:55115/141‑1 (MQ=255)                                                                                                 ttattaCATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGaa                                         >  1:206862/1‑141 (MQ=255)                                                                                                 ttattaCATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGaa                                         >  1:222470/1‑141 (MQ=255)                                                                                                         tttGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATaa                                 >  1:156750/1‑141 (MQ=255)                                                                                                          ttGGAGCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAAt                                <  1:28558/141‑1 (MQ=255)                                                                                                                     tGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGgt                     >  2:108244/1‑141 (MQ=255)                                                                                                                        aaTATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGtt                  >  1:177484/1‑141 (MQ=255)                                                                                                                        aaTATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGtt                  >  2:20955/1‑141 (MQ=255)                                                                                                                           atTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAACTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTggg               <  1:171532/141‑1 (MQ=255)                                                                                                                                   aaaTAGTGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAatgatg       <  1:270772/141‑1 (MQ=255)                                                                                                                                   aaaTAGTGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAatgatg       <  1:202352/141‑1 (MQ=255)                                                                                                                                   aaaTAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAatgatg       >  1:76945/1‑141 (MQ=255)                                                                                                                                      taatGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAATGATGAAc    >  2:89061/1‑141 (MQ=255)                                                                                                                                       aatGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAATGATGAAcc   <  2:142074/141‑1 (MQ=255)                                                                                                                                        atGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAATGATGAACCa  >  2:204551/1‑141 (MQ=255)                                                                                                                                        atGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAATGATGAACCa  >  1:6518/1‑141 (MQ=255)                                                                                                                                        |                                                                                                                                             ATATATTAGACATTTCTATGCAGGTGAAGCTTGAGGGACATGAACAGCCACTTGAATATAAATTAGAAGAAATAGAAACTATGCATTAAGGAGGCATTATTACATTTGGATCGCAATGAAATATTTGAAAAAATAATGCGTTTAGAAATGAATGTCAATCAACTTTCAAAGGAAACTTCAGAATTAAAGGCACTTGCAGTTGAATTAGTAGAAGAAAATGTAGCGCTTCAACTTGAAAATGATAATTTGAAAAAGGTGTTGGGCAATGATGAACCA  >  CP000730/522864‑523139

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