Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,516	C→A	31.5%	Q196K (CAA→AAA)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,516	0	C	A	31.5%	84.9 / 44.0	62	Q196K (CAA→AAA)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (26/16);  new base A (7/12);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 9.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

AAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  USA300TCH1516_ALE/111382‑111650                                                                                                                                        |                                                                                                                                       aaGTTACAGAAGAAAACATAAA‑AAACGCGCTCGAAAACCTCTTGTCACACACACGAGAGAGCAATTTTAAAGATTTAATCAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                                                                   <  1:264173/140‑1 (MQ=9)   gTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                                                                                                 <  1:186073/140‑1 (MQ=255)    tcaaaGAAGAAAACGTAAA‑AAAAGAAATCGAACACTTTAAGTTTTTTGCGCAAGATAGCAATTTTAAAGAGTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                                                                <  2:462002/138‑1 (MQ=18)     tAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATCTAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                                                               >  2:214123/1‑140 (MQ=255)     tAAAGAAAAAAAAAAAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAataagatc                                                                                                                               <  2:550026/140‑6 (MQ=21)      atagcagaaTACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                                                              <  2:381008/135‑1 (MQ=255)       aagaagAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATATCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                                                             <  2:92467/140‑1 (MQ=255)        agaagaAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTaa                                                                                                                            >  2:205004/1‑140 (MQ=255)         gaagaaAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                           <  2:175927/140‑1 (MQ=255)            gaaAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTAATCAGCGTAATATCAATTAACTaa                                                                                                                        >  2:131075/1‑140 (MQ=255)            cgaaaaATAAACAAAAGAAAGCGAAAACTTTAAGTTTTTTGCGCAATAAAGCAATTTTATAGATGTAATGAATTATAAAGGTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                                         <  2:545843/138‑1 (MQ=17)             aaaaactAAA‑AAAAGAAAGCGAAAACTTAAAGTTTTTTGCGCAATATAGCNATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                       <  2:140503/134‑1 (MQ=255)             aaaaCATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                       >  1:401428/1‑140 (MQ=255)                    aaa‑aaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                                                                <  2:242226/140‑1 (MQ=255)                     aa‑aaaaGAAATCGAAAACTCTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATta                                                                                                               >  1:224248/1‑140 (MQ=255)                      a‑aaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                                              >  1:408018/1‑140 (MQ=255)                           aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATg                                                                                                          >  2:433165/1‑140 (MQ=255)                            gAAATTGAAAACTTTAAGTTTTTGGCGCAATATAGCAATTTTAAATATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                                                                         <  2:164254/140‑1 (MQ=255)                              agtCGAGAACTGCAAGGTTTTTGCGCAATATAGTAATTGTAAAGATTTAATGAATTATAAAGATGGGGATATATGATATGATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATTATTATAATGTaa                                                                                                       <  2:173698/138‑1 (MQ=17)                              aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGAAATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                       <  1:147424/140‑1 (MQ=255)                               aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                      <  2:207953/140‑1 (MQ=255)                                gcGAAAACTTTAAGTGTTTTGCGCAATATAGCAATTTTAAAGAATTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                     <  2:69240/139‑1 (MQ=17)                                  gAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                   <  2:368429/140‑1 (MQ=255)                                  gAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                   <  1:80283/140‑1 (MQ=255)                                    aaaCTTTAAGTTTTTTGCGCAATATAGTAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                 <  2:101244/140‑1 (MQ=255)                                     aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                                >  1:484481/1‑140 (MQ=255)                                     aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                                >  2:419426/1‑140 (MQ=255)                                     aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAtc                                                                                                >  1:85838/1‑139 (MQ=255)                                      aCTTGAAGTTTTTAGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                                                               <  2:196022/140‑1 (MQ=255)                                       cTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAGTATAATGTAAAACAATTAc                                                                                              >  2:564909/1‑140 (MQ=255)                                         ttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAACACAATTACgc                                                                                            >  1:63437/1‑139 (MQ=255)                                         ttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgc                                                                                            >  2:168039/1‑139 (MQ=255)                                           aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTCAACCAATTACGTac                                                                                          >  1:123702/1‑139 (MQ=255)                                                tgttGCGGGATAGAGCAATTTTAAAGATGTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGat                                                                                     <  2:547981/138‑1 (MQ=14)                                                      gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATAACa                                                                               >  2:47937/1‑140 (MQ=255)                                                       cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATACAGCGGTGCCGAGTTATTCGGCTAAATATCAATTAACTCAGGATGATTATAATGTAAAACAATTACGCAAAAgatatgatat                                                                              >  2:347031/1‑140 (MQ=14)                                                          tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAcc                                                                           >  1:101335/1‑140 (MQ=14)                                                                aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACATAAAAAATATGACATAACAACaaa                                                                     >  2:273817/1‑140 (MQ=255)                                                                 aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAACGTAAATCAATTACGTAAAAGATATGACATACCAACAAAt                                                                    <  1:297807/140‑1 (MQ=255)                                                                     tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGACAATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAg                                                                >  1:565305/1‑140 (MQ=255)                                                                        aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                             >  1:17825/1‑140 (MQ=255)                                                                        aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                             >  1:318301/1‑140 (MQ=255)                                                                         gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCCACAAATAAAGCGcc                                                            >  1:152633/1‑140 (MQ=255)                                                                               attaaTTACAAAAATGTATAAATATCATATAACCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCt                                                      <  2:524057/137‑1 (MQ=17)                                                                                 gAATTATAAAGATGGAGATATATCATATAATACAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATAAGACATAAAAACAAACAAAGCACCGAAGCagt                                                    >  2:326779/1‑140 (MQ=17)                                                                                               gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCAGTTGTTGAAAGGTACa                                      >  2:266129/1‑140 (MQ=255)                                                                                                agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                     >  1:19536/1‑140 (MQ=255)                                                                                                 gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                    <  1:294527/140‑1 (MQ=255)                                                                                                  atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCAGTTGTTGAAAGGTACAggg                                   >  2:21915/1‑140 (MQ=255)                                                                                                  atatatCATATAAGCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                   <  1:87269/140‑1 (MQ=255)                                                                                                               tCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                      >  1:404279/1‑140 (MQ=255)                                                                                                                 gagagGTGCCGAGTTATTCAGATCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                    <  1:276543/139‑1 (MQ=255)                                                                                                                   gagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                  >  1:179620/1‑140 (MQ=255)                                                                                                                    agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca                 >  1:42363/1‑140 (MQ=255)                                                                                                                    aaacaccGAGTAATTCAGCTCAATATCAATTAACTAATGACGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca                 <  1:80059/135‑1 (MQ=255)                                                                                                                        gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg             >  2:572155/1‑140 (MQ=255)                                                                                                                        accGAGTTAGTCAGCCCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg             <  1:67787/139‑1 (MQ=255)                                                                                                                             gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACACTTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat        >  1:198058/1‑140 (MQ=255)                                                                                                                             gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat        <  2:501657/140‑1 (MQ=255)                                                                                                                             gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGATACAGGGAATGTAAAAGGTTCATCAGTTGGatat        >  2:308058/1‑140 (MQ=255)                                                                                                                              ttATTCAGCTCAATATCAATTAACTAATGATGACTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGAtata       >  2:449546/1‑140 (MQ=255)                                                                                                                                   cAGCTCAATATCAATTAACTAATGATGATTATAATGTAACACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCAACAGTGGGATGTaaaaaa  >  2:300219/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       AAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  USA300TCH1516_ALE/111382‑111650

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